Strongly Inhomogeneous Phases and Non-Fermi Liquid Behavior 
in Randomly Depleted Kondo Lattices 
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We investigate the low-temperature behavior of Kondo lattices upon random depletion of the local 
/-moments, by using strong-coupling arguments and solving SU(iV) saddle-point equations on large 
lattices. For a large range of intermediate doping levels, between the coherent Fermi liquid of the 
dense lattice and the single-impurity Fermi liquid of the dilute limit, we find strongly inhomogeneous 
states that exhibit distinct non-Fermi liquid characteristics. In particular, the interplay of dopant 
disorder and strong interactions leads to rare weakly screened moments which dominate the bulk 
susceptibility. Our results are relevant to compounds like CeajLai.^CoIns and Ce^Lai-^Pba. 



Heavy-fermion metals jj are in the focus of correlated 
electron physics for numerous reasons. They feature a 
strongly correlated heavy Fermi-liquid (HFL) state that 
arises from the Kondo screening of a lattice of localized 
/-moments by conduction electrons. The HFL phase can 
undergo quantum phase transitions to various magnetic 
phases, e.g., spin-density wave and spin-glass states. Re- 
markably, even the Fermi-liquid (FL) regime of heavy- 
fermion metals is not completely understood. Experi- 
mentally, different observables such as entropy, bulk sus- 
ceptibility, and resistivity, are employed to define charac- 
teristic crossover temperatures, but their precise relation 
to microscopic processes is not known. Usually, the be- 
havior cannot be characterized by a single energy scale 
only; this situation is markedly different from that of a 
single magnetic impurity in a metal, where the single- 
impurity Kondo temperature is the only low-energy 
scale, and all thermodynamic observables are universal 
functions of TjT^ only. 

A fascinating modification of Kondo lattices is the de- 
pletion of the local-moment sector, by randomly replac- 
ing magnetic by non-magnetic rare-earth ions. This al- 
lows one to continuously tune the system between the 
dense Kondo lattice and the isolated impurity - both 
limits are known to exhibit FL behavior (provided that 
no magnetic or superconducting order intervenes). This 
raises fundamental questions: (i) Do the two FL phases 
evolve continuously into each other upon changing the 
concentration of /-moments? (ii) Over what range of 
doping is a picture of independent impurities applicable? 

Early experimental studies on Ce^Lai-^Pbs indi- 
cated that a single-impurity picture may be valid over 
a rather large range of /-moment concentration: it was 
found that even at x = 0.6 both transport and ther- 
modynamic measurements could be described using a 
single-impurity Kondo model down to the lowest temper- 
atures studied (one-tenth of the estimated Tjp). These 
results are particularly surprising because the "Kondo 
screening clouds" should be strongly overlapping in this 
regime. In contrast, recent detailed investigations of 
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FIG. 1: Schematic phase diagram of the depleted Kondo 
lattice. In the dense limit, n/ < 1, a coherent heavy FL is 
formed below T co h (dashed), whereas in the dilute limit, n/ — * 
0, single-impurity physics prevails, with local FL behavior 
below (dash-dot). There is a large regime of intermediate 
n/ where low-energy properties are strongly inhomogeneous, 
and FL behavior is suppressed, a) For weak inter-moment 
exchange we have T co h < \ whereas b) stronger inter- 
moment exchange enhances T co h and shifts the boundary of 
the inhomogeneous regime towards smaller n/. 



Ce^Lai-ajCoIns Q interestingly revealed different prop- 
erties: FL behavior, as inferred from the resistivity, is 
quickly suppressed to very low T when moving away both 
from the dense and dilute limit. The energy scale char- 
acterizing the onset of moment screening in the suscep- 
tibility is enhanced by more than an order of magnitude 
in the dense case compared to the dilute limit. 

On the theoretical side, the behavior of depleted Kondo 
lattices has been studied rather little. The experiments 
on Ce^Lai-^CoIns triggered an interpretation in terms 
of a phenomenological two-fluid model Q (see Ref. for 
a related mean-field approach). However, detailed mi- 
croscopic investigations taking into account the dopant 
disorder (apart from averaged treatments using CPA @) 
arc lacking. 

The aim of this Letter is to close this gap: using 
both strong-coupling arguments and detailed numerical 
results, obtained from solving SU(7V) saddle-point equa- 
tions on finite lattices, we shall present a consistent pic- 
ture of the evolution of the HFL upon random /-moment 
depletion 0,13 . 
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Our main findings can be summarized as follows. The 
crossover from HFL to the single-impurity FL takes place 
at a concentration of / moments that strongly depends 
on the number of delocalized electrons per unit cell, n c . 
In particular, at large n c single-impurity behavior is ob- 
served even for dense / moment concentrations (as in 
Ce^Lai-^Pbs). Between the two FLs we find a wide 
crossover regime with strongly inhomogeneous states at 
low temperature: The local density of states (LDOS) in 
the conduction band near the Fermi level shows signa- 
tures of localization, associated with island formation, 
and the uniform susceptibility does not saturate down 
to very low temperatures due to the presence of poorly 
screened moments. Importantly, the broad distributions 
of local quantities, such as the local susceptibility xioc, 
and the associated non-FL behavior do not originate from 
disorder in the bare LDOS or in the local Kondo cou- 
plings (as in so-called Kondo disorder models [§]), but 
arise through the collective screening of the randomly 
positioned /-moments by the conduction electrons. Our 
results are summarized in the phase diagrams, Fig. H] 

Model. The standard Kondo lattice model with ran- 
domly depleted /-moments reads 
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where c ka represent the conduction electrons with mo- 
mentum fc, spin a, a band width D, and a band filling 
of n c (n c = 1 corresponds to half filling, in the follow- 
ing we shall assume n c < 1). The S r are the spin-1/2 
moments on the Nf occupied sites {r} of the / electron 
lattice, with a concentration of nf = Nf/N s where N s is 
the number of conduction electron sites [rif = 1 is the 
dense Kondo lattice). 

Strong- coupling analysis. In the limit Jk 3> t in 
Eq. (JJJ the c-electrons and /-spins form local singlets. In 
the dilute limit, n/ -C 1, one ends up with Nf local sin- 
glets which act as potential scatterers for the remaining 
c-electrons - the system is a FL In the dense limit, 

rif — 1, it is useful to think in terms of holes: after form- 
ing N s singlets we are left with N s (l — n c ) mobile holes 
in the conduction band which are HFL quasi-particles: 
clearly this limit is also a FL. Now, as the /-moments 
are randomly depleted, the band of HFL quasi-particles 
is also depleted. Let us focus on nf = n c : In the strong- 
coupling limit all conduction electrons are bound into 
localized singlets, and the system has a gap of order Jk 
to all excitations - this is apparently not a Fermi liquid. 
This discussion implies that the ground state of Eq. Q 
evolves discontinuously upon tuning of nf. In particular, 
n f = n c appears to be the dividing line between HFL and 
single-impurity FL behavior. In the region rif « n c , we 
expect that the system is most sensitive to randomness 
(introduced by /-spin depletion), such that localization 
phenomena and strong inhomogeneities are likely to oc- 
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FIG. 2: (color online) Summary of the mean-field results: 
(a) T* — (br)/D as function of the local-moment concentra- 
tion n/ for Jk I D = 1.25 and different band fillings n c . Note 
the sharp change around n c = rif. The inset shows the distri- 
bution of br/(br) for the different n c at fixed nf =0.5. The 
strongly inhomogeneous regime is characterized by a bimodal 
distribution, here most pronounced at n c = 0.4. (b) Evolution 
(from top to bottom) of the mean-field density of states from 
the HFL to the dilute impurity limit at fixed n c = 0.4. 



cur. These general arguments will be supported by our 
mean-field calculations below. We note that the gap at 
n f = n c (of similar origin as the gap of the Kondo insu- 
lator, nf = n c = 1) is a feature of the strong-coupling 
limit, and not expected to survive for Jk t. 

Mean-field theory. To obtain quantitative results, we 
now turn to a solvable large- N limit of Eq. |JTJ that al- 
lows us to access arbitrary JK/t. Using slave fermions 
to represent the local moments, S r = \f}. a <r aa > f ra i with 
frafra = 1> decoupling the Kondo interactions with aux- 
iliary fields b r , and neglecting temporal fluctuations of 
the b ri we obtain the mean-field Hamiltonian 
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Note that H m { represents the N = oo saddle-point so- 
lution of the SU(7V) Kondo lattice model, with a fully 
antisymmetric representation of the local moments. The 
mean- field parameters b ri ^i r are determined by the self- 
consistency conditions 



i = (£/™>. 2^ = ^(4,/™). 



(3) 



At high temperatures the only solution is b r = fi r = 0, 
corresponding to quasi-free local moments with Curie 
susceptibility. Upon lowering T, all b r will become non- 
zero at Tft \ the (mean-field) single-impurity Kondo tem- 
perature, which is determined by Jk and the LDOS of 
the conduction band only. Hence, the onset of Kondo 
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screening at 01 does not depend onn/ and on the 
particular disorder realization. Within mean-field theory, 
the simplest energy scale describing the low-temperature 
behavior is T* = (\b r (T — 0)\ 2 ) / D, where (...) denotes a 
disorder average. For a single impurity, the ratio T* /T& 
is a constant of order unity. In contrast, in the dense 
Kondo lattice T*/T^ < 1 depends on the band filling 
n c (but not on Jx at weak coupling); T* can be associ- 
ated with the coherence temperature T cn h below which 
the Kondo lattice behaves as a FL hi, 121. 

Numerical results. We have performed large-scale nu- 
merical simulations of the mean-field equations on square 
lattices of up to 20 x 20 sites for different n c and n f . 
Fig. |2Ia) shows a summary of the mean-field condensa- 
tion parameters. Upon increasing rif the scale T* stays 
roughly constant until nj « n c and then it drops sharply 
towards its nf — 1 value. This supports the picture that 
emerged from strong-coupling: rif ~ n c separates HFL 
from single-impurity FL behavior. The distribution func- 
tion of (b%) becomes bimodal when nr « n c pointing to 
strong inhomogeneities at this filling. Fig. H£b) shows 
the evolution of the mean-field density of states of the 
HFL upon random /-moment dilution. For the dense 
lattice the Fermi level lies in region of the lower band 
with a large density of states. The number of states in 
this heavy-electron band is reduced upon /-moment de- 
pletion. Near nt = n c the two bands merge, and the 
Fermi level lies near the edge of the bands where the 
states are most localized. (We have studied the inverse 
participation ratio of the mean-field wave functions [l^j 
and confirmed this expectation.) The system at this fill- 
ing is thus in a strongly inhomogeneous phase. 

To gain a better understanding of the underlying mi- 
croscopies, let us focus on local properties of a partic- 
ular disorder realization. Fig. [3] shows the spatial dis- 
tributions of the local magnetic susceptibility xioc and 
the LDOS A(u) (defined as the imaginary part of the 
local c electron propagator) at the Fermi level (aver- 
aged over a small energy interval to avoid discretiza- 
tion effects). Let us first discuss the limiting cases: in 
the dilute limit, < 1, we expect a single- impurity 
picture to hold: Each moment is individually Kondo- 
screened, and the LDOS is suppressed at the f -moment 
locations because of the resonant scattering. For the op- 
posite limit, 1 — rif -C 1, we can develop a single- vacancy 
picture: the missing / moments act as scattering cen- 
ters in the effective two-band system, and the LDOS is 
now suppressed at the vacancy locations. This picture 
is realized in Fig. [21 with n c — 0.8 fixed and varying nj. 
Surprisingly, a single-impurity behavior (with negative 
correlation between LDOS and the site occupancy of /- 
moments) is realized even at large /-moment concentra- 
tions like rif = 0.5. It is only at rif = 0.95 that one finally 
sees the single-vacancy behavior (positive correlation be- 
tween LDOS and site occupancy); this is consistent with 
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FIG. 3: (color online) Spatial distribution of the slave boson 
amplitudes (b r , left), the low-energy LDOS (A r (u — 0), mid- 
dle), and the local susceptibilities (xioc, r(T= 10 -3 ), right) for 
different n f = 0.05, 0.5, 0.95 at a band filling of n c = 0.8. The 
Kondo coupling is Jk/D = 1.25. 



rif w n c separating the HFL and single- impurity FL. 
Fig. E] illustrates that varying n c can also be used to 
tune between the single-impurity and single- vacancy lim- 
its (here for a fixed realization of disorder!). 

We now turn to the strongly inhomogeneous phases 
occurring for intermediate values of rif around rif = n c . 
Fig. EJs clearly shows island formation in the low-energy 
LDOS. This in turn gives rise to huge fluctuations in 
the Kondo temperatures of individual moments, result- 
ing in broad distributions of, e.g., xi oc and the existence 
of weakly screened moments down to lowest T . Similar 

a)LDOS,n t . = 0.1 b) LDOS, n c = 0.2 c) LDOS, n c = 0.8 




FIG. 4: (color online) Spatial distribution of the LDOS, 
A r (u> = 0), for a specific disorder realization with Uf = 0.2, 
but for varying conduction electron filling n c = 0.1,0.2,0.8. 
The crossover from the single-vacancy (a) to the single- 
impurity (c) regime can be nicely seen, with a tendency to- 
wards island formation in between. The /-moment locations 
can be identified from the dark sites in (c). 
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FIG. 5: (color online) Temperature dependence of x an d 
mean-field parameters, (a) /-moment susceptibility Xff> m 
the inhomogeneous regime nf ~ n c — 0.4, \ff does not satu- 
rate to the lowest T studied, (b) Local susceptibility, xioc(T), 
at three different sites, illustrating the broad range of behav- 
ior. (c,d) Mean-field parameters b r (T) and j-i r (T) at the same 
sites. Jk/D — 1.25 and for (b,c,d) nf — n c = 0.4. 

to impurities in disordered metals and disordered Kondo 
lattices , these broad distributions of Tk produce non- 
FL behavior in both thermodynamic and transport prop- 
erties - these arguments qualitatively apply to our situ- 
ation as well. 

In Fig. [5] we show numerical data illustrating the be- 
havior of the susceptibility; other observables (entropy, 
specific heat etc.) will be presented elsewhere [l4j . 
Fig. Eta) shows how Xff(T) in the Kondo lattice and 
impurity limits saturate at low-T. In contrast, for inter- 
mediate nj the Xff{T) does not saturate to the lowest 
T studied (finite-size effects prevent calculations at lower 
T). While we cannot rule out saturation of Xff at ex- 
tremely low temperatures, it is clear that FL behavior is 
strongly suppressed in this regime, as illustrated in the 
phase diagram, Fig. [J- this behavior is very reminiscent 
of the experimental data on Ce x Lai_ x CoIn5 Fig.|S{b) 
shows Xloc(T) for three different sites, demonstrating 
that xioc is strongly inhomogeneous (while the b r , \x r are 
more homogeneous), and the behavior of x//(T) arises 
from a fraction of the / moments only. 

Inter-moment correlations. So far, we have neglected 
the effects of a magnetic inter-moment interaction, I, ei- 
ther due to direct exchange or mediated by the conduc- 
tion electrons. These interactions can be important for 
not too small nf ||, provided that their strength is com- 
parable to or larger than Tk- (Note that / ;§> Tk usually 
leads to a magnetically ordered state.) In the dense limit, 
the coherence temperature will be enhanced compared to 
the I — situation, as the local moments are quenched 
by / due to inter-moment correlations. To assess the 
effect of I for intermediate nj we have performed cal- 
culations using the large- N formalism of Ref. 0. Our 
preliminary results (not shown) indicate that the rare 
unqucnchcd moments tend to be suppressed for larger 
n/ and sufficiently strong /, but islands and broad Xioc 
distributions survive for smaller n / over a sizable param- 



eter range . The resulting phase diagram is sketched 
in Fig. 03, and is qualitatively remarkably similar to the 
one of Ce x Lai_ x CoIn 5 0. 

Conclusions. We have studied randomly depleted 
Kondo lattices over the full crossover from the dense lat- 
tice to the dilute impurity limit. We found that n/ = n c 
appears to be the dividing line between the two limit- 
ing FL behaviors; in the vicinity of rif = n c the sys- 
tem is generically strongly inhomogeneous and displays 
non-FL behavior. Our results are relevant to a vari- 
ety of Kondo-lattice metals, e.g. Ce K Lai_a;CoIn5, that 
show both HFL and single-impurity FL behavior as the 
/-moments are randomly depleted, without signs of mag- 
netic order. While the condition n/ ~ n c likely only ap- 
plies to heavy-fermion metals described approximately 
by the 5=1/2 single-channel model, Eq. 0J, we ex- 
pect the presence of two FL regimes separated by inho- 
mogeneous phases (Fig. ^| to be a generic feature of /- 
moment-depleted heavy FLs. We speculate that a rather 
large (effective) n c is the reason why, e.g., Ce x Lai_ a Pb3 
displays single-impurity like behavior for sizeable range 
of x = nf . 
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